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Abstract. We discuss the well known Einstein and the Kubo Fluctuation Dissipation 
Relations (FDRs) in the wider framework of a generalized FDR for systems with 
a stationary probability distribution. A multi-variate linear Langevin model, which 
includes dynamics with memory, is used as a treatable example to show how the 
usual relations are recovered only in particular cases. This study brings to the fore 
the ambiguities of a check of the FDR done without knowing the significant degrees 
of freedom and their coupling. An analogous scenario emerges in the dynamics of 
diluted shaken granular media. There, the correlation between position and velocity of 
particles, due to spatial inhomogeneities, induces violation of usual FDRs. The search 
for the appropriate correlation function which could restore the FDR, can be more 
insightful than a definition of "non-equilibrium" or "effective temperatures" . 
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1. Introduction 

The idea of a link between dissipation and fluctuations dates back to Einstein with 
his work on the Brownian motion [1] and his relation between mobility (which is a 
non-equilibrium quantity) and diffusion coefficient (which is an equilibrium quantity). 
Later Onsager [2j 13] with the regression hypothesis and Kubo [HE], with linear response 
theory, investigated in a deep way the issue of the fluctuation-dissipation relation (FDR). 
In the last decades there has been a renewed interest in this topic, see the contributions 
of Morris, Evans, Cohen, Gallavotti and Jarzynski [61 [3, [8] to cite just some of the most 
well known attempts (for a recent review, see [9]). 

The FDR theory was originally introduced in the Hamiltonian systems near 
thermodynamic equilibrium. However it is now clear that a generalized FDR exists, 
under very general assumptions, for a large class of systems with a "good statistical 
behaviour" , i.e. with a relaxation to an invariant (smooth) probability distribution. We 
stress that such a condition is quite common, e.g. in any system with a finite number 
of degrees of freedom whose evolution rules include some randomness (for instance non- 
linear Langevin equations). Unfortunately the explicit form of generalized FDR depends 
on the shape of the invariant probability distribution (which is typically unknown), 
however this is only a technical difficulty without conceptual consequences [9]. 

We are not concerned with systems without a stationary probability distribution, 
e.g. systems showing aging [TO], E] and glassy behaviour [12J. On the contrary, we 
consider here systems with an invariant probability distribution satisfying the hypothesis 
for the validity of the generalized FDR mentioned before. For such systems, in our 
opinion, there is some confusion about what form of FDR has to be expected. When 
couplings between the chosen observable and other degrees of freedom are ignored, the 
wrong FDR is expected, i.e. the response function is compared to the wrong correlation: 
this leads to what is often called a "violation" of FDR. 

The structure of the paper is the following: at first, in Section 2, we discuss the 
different kinds of FDR obtained in the statistical mechanics framework. In Section 3 
we analyze a one- dimensional Langevin equation with memory, which can be mapped 
to a multivariate Langevin equation without memory: this example can be worked out 
analytically and well illustrates all the main points of our discussion. In section 4 some 
numerical results on a driven granular gas model are discussed to support our general 
discussion with physical examples. Section 5 is devoted to concluding remarks. 

2. Linear response in statistical mechanics 

Let us briefly recall three different kinds of Fluctuation Dissipation Relations (FDR), 
commonly used in statistical mechanics when a small impulsive perturbation is applied 
to a stationary system. These three formulae share the feature of relating the system's 
linear response to an appropriate two-time correlation computed in the unperturbed 
system. Anyway these relations have different fields of application and must be adapted 
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with care. For a more pedagogical introduction, we first set the notation, and then we 
discuss the three FDR versions. 

2.1. Linear response functions 

We adopt the following notation: 8(t) denotes the Dirac delta function, and 8ij 
the Kronecker delta, we use the overline ~ for non- stationary averages over many 
realizations and (•) for averages using the unperturbed stationary probability in phase- 
space (assuming ergodicity, this is equivalent to a time-average over a long trajectory). 
Accordingly, we use the shorthand notation 

Cab® = (A(t)B(0)) (1) 

to denote the two-time correlation function between observables A(K(t)) and B(X(t)), 
with X(t) the state of the system at time t. Let us introduce the matrix of linear 
response functions, whose ij element reads 



i.e. the mean response of the variable X{ at time t to an impulsive perturbation applied 
to a variable Xj at time 0. If the dynamics of the system is given, for instance in the 
form ^ = f(X), the mean linear response of the 2-th degree of freedom to a small 
perturbation of the j-th component of the vector field fj — > fj + 5fj can be expressed 
as 



8X i (t) = dsR XiXj (t-s)8f J (s). (3) 

J — oo 

The case of an impulsive perturbation at time 0, 5fj = 8Xj(0)8(t), gives back the 
definition 

Let us also consider the historically important case of a Hamiltonian system at thermal 
equilibrium: in this case the perturbation is typically defined on the Hamiltonian, i.e. 
7i — > 7i + 8H, with 8H = —8h(t)B(K). A linear response function of an observable A 
to an impulsive field 5h(0)8(t) at time t — 0, R.Ah(f) = f^y determines the behaviour 
of 8A(t) for a generic small perturbation 8h(t): 

~8A(t)= I dsR Ah (t-s)8h{s). (4) 

J — oo 

There are not conceptual differences between the two procedures to perturbe the state 
(i.e. with a 8Xi(0) or the introduction of an extra term in the Hamiltonian): we can 
simply consider A and B as two variables of the system. For instance, consider the case 
where X is the position of a colloidal particle evolving according to 




with 7]i independent normalized white noises, i.e. Gaussian processes with (t^) = 
and (ili(t)r]j{t')) = 8ij8(t — t'). Here, the effect of a perturbation AH = —Xj8h(t) is 
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equivalent to a 5 fa = Sij^^- and therefore, from a comparison of Eqs. ([3]) and (jl]), one 

has R Xl h = Rxixjl- 

2.2. Three different Fluctuation Dissipation Relations 
We can now discuss the three forms of FDR we are interested in: 
(i) the generalized FDR, denoted GFR in the following 

flW*) = c W*). ( 6 ) 

where Sj, see below, depends on the invariant distribution density in phase space. 
This relation is valid (under quite general assumptions, see [13], El [El [161 EH E]) 
in a dynamical system, whose state is completely determined by the phase space 
coordinate X, and whose dynamics induces an invariant measure of phase-space 
p(X) | and 

dlogp(X) 
so that, the correlation reads: 



dx j 

(ii) the so-called Einstein relation, in the following referred as EFr[§|, 

R (t) ~ M{t) - ° AA{t) 

The most known example of the above relation is given by a Brownian colloidal 
particle diffusing in an equilibrium fluid at temperature T. Such a system is 
described by a linear Langevin equation 



mv = —'-/v 



+ V^fTv (10) 



where rj is a normalized white noise, i.e. a Gaussian process with (77) = and 
(r](t)r](t')) = 5(t — t'). It is easy to show that after a small impulsive force 
5h(0)5(t) at time 0, the velocity perturbation, which at t = is 5v(0) = 5h(0)/m, 
decays as 5v(t) = (k>(0)e~ 7 *. On the other hand, since C vv {t) = C„t,(0)e~ 7 *, where 
C vv (0) = (v 2 ) = T/m, the EFR © holds, with A = v (see [TBI). The Green- 
Kubo relations, relating transport coefficients to the time-integral of unperturbed 
current-current correlations are the extension of relation (Q to generic transport 
processes; 

I more precisely one assumes that p(X) is a smooth non-vanishing function. This condition surely 
holds if some noise is included in the dynamics. 

§ for simplicity we adopt the name "Einstein relation" which, in the literature, has been typically used 
to denote the time-integral of relation (|9|), e.g. the formula fi = [3D relating the mobility /i to the 
self-diffusion coefficient D. 
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(iii) the classical Kubo relation, hereafter denoted as KFR, 
Rah = —f3-j-CAB(t)- 



(11) 



This relation holds - for instance - when a Hamiltonian system, whose statistics 
is described by the canonical ensemble with temperature T = 1//3, is perturbed 
by a Hamiltonian variation 57i = —5h(t)B, which defines the perturbing force 
h and the corresponding conjugate field B (see [5], [9]). Relation (TTTT) also holds 
for Langevin equations with a gradient structure, e.g. Eq. (jHJ). Introducing the 
quantity Xab(£) = C RAh(t')dt' [jj], from (fill) one has 



In some literature pj2,[20j[2T] a deviation from (fT2l) is indicated as a failure of FDR, 
a mark of being far from equilibrium, and is used to define new "non-equilibrium" 
temperatures. 

As the name suggests, the GFR (iii) includes both EFR (i) and KFR (ii) for some 
choices of the dynamics or of the stationary distribution p(X) [9J. This can be shown 
for the case A = Xj and B = Xj and is easily generalized to any other discussed 
above: 

• GFR — > EFR: this happens when the invariant distribution in phase space is of the 
form p(X) = exp(— ^ X 4 2 )/Z (being Z a normalizing constant), then one has 
Sj = Xj/T and it is immediate to obtain EFR starting from relation GFR; 

• GFR — > KFR: this happens for Hamiltonian systems in the canonical ensemble, 
or Langevin equations with gradient structure when a small force is applied: in 
both cases one has p(X) = exp[— /37i(X)]/Z, so that the GFR involves the quantity 
Sj = P§^~- If the dynamics is given, for instance, by an overdamped Langevin 
equation of the kind ([5]), one has that Sj = —(3^Xj, i.e. RxiX 3 = Pl{Xi(t)Xj(Q)) 
and, considering the discussion after Eq. ([5]), the KFR Rxih = P(Xi(t)Xj(0)) is 
immediately derived. 

3. Langevin equation with memory 

Let us now consider a system that does not necessarily fall neither under the hypothesis 
of EFR, neither under those of KFR. Our choice here goes to a linear stochastic equation 
with memory: 



where m is the mass of the tracer, k is the constant of an elastic force, F(t — t') is the 
friction kernel, r] (t) is a stochastic force acting as a thermostat. 

|| For this last quantity we use a notation where the second subscript B directly refers to the observable 
conjugate to the perturbed field h: even if it is not self-evident, this has the advantage of being coherent 
with the notation widely used in the literature. 



XAB(t)=PC AB (0)-(3C AB (t). 



(12) 




(13) 
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When there is no memory, i.e. T(t — t') = 2 r y5{t — t') and (i](t)r](t')) = 2 , ~fT5(t — t'), 
the usual Langevin equation is recovered, leading to EFR in the case k — 0, Eq. ffTUj) . 
or KFR in the overdamped case, Eq. (jSJ). 

For a generic memory kernel, Kubo [5] has also shown that EFR or KFR are 
recovered, provided that 

(r ) (t)r ] (t'))=TT(t-t'). (14) 

Here we consider a memory kernel of the kind: 

r(t-t) = 2 7/ 5(t-t / ) + -exp f-^n (15) 

and a noise ^(t) = p/(t) + p s (£) where p/(£) and p s (t) are two independent Gaussian 
processes with zero means and 

( Pf (t)p f (t>)) = 2T fV 5(t - t') ( Ps (t) Ps (t')) = T s ^e-^. (16) 

r s 

so that condition <HM is recovered when Tf = T s = T ("s" and "f" subscripts stay for 
"slow" and "fast" respectively). In general, however, Tf ^ T s . Note also that Tf and 
T s have the dimension of a temperature: indeed the model in its overdamped limit has 
been proposed as an example of system coupled to two different baths acting on different 
time-scales H91I201. 



3.1. A Markovian equivalent model 

The first observation about the system in study is that, because of the memory term, 
its dynamics after time t cannot be deduced by the knowledge of x and v at time t: in 
fact, the evolution depends on its history, i.e. the dynamics is non-Markovian, and the 
GFR cannot be directly applied. 

However, it is possible to recover Markovianity, at the price of adding additional 
degrees of freedom. In other words, it is possible (and it will be done in the next section) 
to recast equation ( fl3|) in a linear, multi-dimensional, Langevin equation: 

<ix , . 

where x e <ft are N-dimensional vectors and A is a real N x N matrix, in general not 
symmetric. In addition, now <f>(t) is a Gaussian process, with covariance matrix: 

(<f )t (t')<f )l (t)) = 5(t-t')D lv (18) 

and the real parts of A's eigenvalues are positive. The stationary probability density 
is [221: 



P(x) = (27r)-^[Deta]-^exp j—^^^J , (19) 

where a is a symmetric matrix determined by the following relation: 

D = Ao + oA T . (20) 
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We can now explicitly study the fluctuation and response properties of the system 
since the dynamics, being now Markovian, satisfies the hypothesis of applicability of 
GFR. First, we recall the definition of correlation matrix Cy(t) = (xi(t)xj(0)), in the 
stationary state, which is time-translational invariant. Then, using the equation of 
motion, it is immediate to verify that C(t) = —AC(t), with initial condition given 
by the covariance matrix between degrees of freedom at equal time: C(0) = a. The 
corresponding solution is: 

C(t) = exp(-At)a (21) 

Note that, in general a and A do not commute. It is also straightforward to recover the 
response function R(t) = exp(— At), since the GFR imposes the following equation to 
hold: 

R(t) = C(t)o- 1 . (22) 
In general this function can be written as: 

^i? Q exp(-A Q t) (23) 

a 

where R a are constant matrices, and A Q are the eigenvalues of A. The i,j element of the 
matrix R(t) is the response function Rxi,Xj(t) for the corresponding degrees of freedom. 
However, at odds with EFR and KFR, this quantity cannot be expressed in terms of 
the correlation Cx it Xj(t) only, since in general all the degrees of freedom are coupled: 

^x,(t) = ]T(a-%CW^ (24) 

k 

which appears as a violation of EFR or KFR, even if GFR is still valid. 

In the following, we will explicitly walk through this analysis in two different limit 
conditions: 

(i) the free case, when the harmonic force kx can be neglected; 

(ii) the overdamped case, when inertia mi can be neglected. 

3.2. Dynamics of the free particle 

In the limit k — 0, and setting m = 1 without loss of generality, equation ( [131) becomes 
v = - lfV -ll f e -^fv{f)dlf + p f {t)+p a (t), (25) 

Ts J—oo 

where we have introduced velocity v = dx/dt. 
Eq. fl25l) can be mapped to ffTTj) by: 

v \ j_ ( y/ 2 7f T fM t) 
^v^X0 2 (t) 





u 



(26) 



where <j)i(t) and fait) are independent normalized white noises, and u(t) is an auxiliary 
variable: 



u(t) = - f e~*£ lv(t') + J^Mt')) ) dt'. (27) 
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Denoting a = Jf/js, v = t s j s , ( 1 = (1 + a)(l + va) and AT = T s — Tf, a 
straightforward calculation of covariance matrix gives: 

° ~ { -a(AT T s /v - a(AT ) ' (28) 

The more general formula for the response as a function of correlations is given by the 
GFR, Eq. (JED: 



= ^-i 1 (v(t)v(0)) + a£ (v(t)u(0)) . (29) 

We can observe two different scenarios: in the case Tf = T s = T, i.e. AT = 0, a is 
diagonal with <j\\ = T and 022 = T/v, condition (fl4l) is restored and, independently by 
the values of other parameters, a direct proportionality between C vv and R vv is obtained 
(EFR). This is not the only case for this to happen: e.g. for fixed v = r s 7 s , one has two 
possible limits: 

'-(? r >)«'-(o' 

respectively for 7^ 7 a , i.e. a — * 0, or 7^ ^> 7 S , i.e. a — > 00. More in general, when 
Tf 7^ T s , the coupling term a 12 differs from zero and a "violation" of EFR emerges 
between the coupling of different degrees of freedom. 

The situation is clarified by Fig. [H where EFR is violated and the GFR holds: 
response R vv (t), when plotted against C vv (t), shows a non-linear and non-monotonic 
graph. Anyway a simple linear plot is restored when the response is plotted against the 
linear combination of correlations indicated by formula (j29|) . In this case it is evident 
that the "violation" cannot be interpreted by means of any effective temperature: on 
the contrary it is a consequence of having "missed" the coupling between variables v 
and u, which gives an additive contribution to the response of v. 

It is interesting to note that, even when Tf = T s (i.e. o\i = cr 2 i — 0), 
C vv {t) oc R vv {t) is a linear combination of two different exponentials, see Eq. (123]) . 
then its derivative is in general not proportional to R vv , i.e. the KFR does not hold. 
An example of this situation will be discussed in Section H] and Figure [H 

The above consideration can be easily generalized to the case where several slow 
thermostats are present: let us suppose that there are N — 1 thermostats at temperature 
T (the fast one and N—2 slow ones) and one at temperature T\. In this case it is possible 
to show that the off-diagonal terms in a are proportional to (T — T x ). 

It is useful to stress the role of Markovianity, which is relevant for a correct 
prediction of the response. In fact the marginal probability distribution of velocity P m (v) 
can be computed straightforward from (|25|) and has always a Gaussian shape. By that, 
one could be tempted to conclude, inserting P m (v) inside GFR, that proportionality 
between response and correlation holds also if Tf 7^ T s , in contradiction with ff28l) . This 
conclusion is wrong, as stated at the beginning of this section, because the process is 
Markovian only if both the variable v and the "hidden" variable u are considered. 
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Figure 1. Free particle with viscosity and memory, whose dynamics is given by 
Eq. (|25[) : here we show the parametric plot of the velocity response to an impulsive 
perturbation at time 0, versus two different correlations. The Einstein relation, which 
is not satisfied, would correspond to a linear shape with slope 1 for the red dashed 
line. The black line shows that the GFR holds. 



3.3. Overdamped limit with harmonic potential 

When k ^ 0, in the overdamped limit we can neglect the left-hand side term in (fl3|) . 
and the equation, after an integration by parts, reads: 

jfX=-(k + —)x+-^[ e^ x(t')dt' + Pf (t) + Ps {t). (31) 



This model has been discussed in the context of driven glassy systems |19j [20] . 
In order to restore Markovianity, we can map this equation in: 




A. _|_ 7« 7s 

7/ IfTs JfT £ 

_J_ 1_ 

T 3 T a 




(32) 



where, as before, 0i and 02 are independent white noises, while the auxiliary variable u 
now reads: 



u(t) = 1 ^ e-^r L(t') + r s ^Mt'))J dt'. (33) 

Equation (1321 describes the overdamped dynamics of the model depicted in Fig. [2j 
a first particle at position x is coupled to a thermostat with viscosity 7/ and temperature 
Tf, and to the origin by a spring of elastic constant k; a second particle at position u 
is coupled to a thermostat with viscosity 7^ and temperature T s ; the coupling between 
the two particles is a spring of elastic constant k! = 7 s /r s . The first particle, when 
uncoupled from the second particle, has a characteristic time 7/ —Jf/k. 
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Figure 2. Sketch of a simple model which is described by Equation (|32|) . 

Using non-dimensional parameters a = 7//7 s and 77 = t s /tj, the covariance matrix 
<t reads 

(Tl _i_ J_AI Ik _ iat \ 
t or t k V k I fQ/A 

Zk _ 1 AT 1 1 AT I • \°^) 

k r k \ k') k r k J 

where r = 1 + - + 77 and AT = T s — Tf. As we can see, a diagonal form for the matrix 



a is not recovered, even for AT = 0, i.e. for this model the EFR never holds. 

Let us consider now the KFR, with a force field h coupled to the variable x. First, 
we note that, from ([22]), 



C = -RAa (35) 
(where commutativity between R and A has been used). Therefore, in general, since 
§j| = Rxx/lf, one has 

M*) _ ( A °)xx ^ ( A °)ux ^ (1f\\ 

dn{0) 7/ 7/ 

Then, it is easy to see that the condition to have KFR is (Aa) ux = 0. Here the matrix 
Aa reads: 

/ ±Tj_ ^at \ 

Aa = 2-k 7^ , (37) 

\ r/r fc T S k' / 

and the condition (Acr) ua; = is equivalent to AT = 0. 

As observed in Fig. [3j in analogy with the previous Fig. (TJ the GFR always holds 
(solid lines), while the KFR is not verified because it ignores the coupling between 
relevant degrees of freedom. In general the local slope —s(t) of the parametric curve 

Xxx 

(t) vs. C xx (t) is given by 

dC xx (t) 



s fA = dXxx = dx ^ 



(3f 



dt 

which, for the model in Eq. (|32[) . has two different limits: 

l/s{t) -> T f for t -> 0, (39) 

l/s(t) -> 1/soo = T s - KAT for i -> 00, (40) 




Correlation 

Figure 3. Overdamped motion of a particle with harmonic potential, viscosity and 
memory, whose dynamics is described by Eq. (|32[) : we show here the parametric plot 
of integrated response versus two different self-correlations. The Kubo formula KFR, 
which is not satisfied, would correspond to a straigth line of slope —1 for the blue 
dashed curve. The black line shows that the GFR holds. 



where 

2r v ; 

Actually, as observed in [El[20], when r s ^> r/ the parametric plot Xxxif) vs. C xx (t) 
takes the form of a broken line with two slopes: the point where the slope abruptly 
changes, corresponds to the intermediate plateau of C xx {t) (i.e. when r s t 3> r/) and 
is located at a position on the x _ax i s ~ 2/oX(°°) with y = : it becomes visible in 
the plot if 2/0 ~ | (i.e. it is not close to 1 or 0). In this case, since K — > 0, the two 
slopes are 1/7/ and 1/T S . This observation has driven a series of real and numerical 
experiments where the parametric plot Xxxif) vs. C xx (t) (or their Fourier transforms for 
the frequency-dependent susceptibilities) were measured for some degree of freedom in 
slowly driven [23] or aging [2H EE] glassy systems, including models for densely packed 
granular materials [26j. In the next section we discuss this limit and other interesting 
cases from the point of view of the GFR: this can be useful to understand when the 
KFR-inspired parametric plot is meaningful and why. 

3.4- Phenomenology of the x vs. C parametric plot 

In this section we probe the hypothesis that the relative relevance of the contributions 
(Aa) ux an( j CAg)*g C xx (t) to the response depends on the time-scale of 

observation. In particular we consider the time-integrals of these two contributions, 
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20 


40 


30 


20 


2 


2/3 


47.3 


12.7 


b 


2 


0.6 


200 


1 


200 


1 


1 


1 


400 


0.5 


c 


2 


0.6 


100 


100 


2 


1000 


0.1 


50 


2000 


1 


d 


10 


2 


1 


50 


10 


50 


1 


0.1 


51.2 


9.76 



Table 1. Table of parameters for the 4 cases presented in Figures 0] and [5j The 
effective time of the "fast" bath is defined as Tf = "Jf/k, while the effective spring 
constant coupling x with u is defined as k 1 — 7 s /r s . 



such that Xxxif) = Qxx(t) + Qxu(t), see Eq. (J36J) : 

Qxx(t) = [0^(0) - C xx (t)} (42) 

7/ 

Q,n(t) = ^^40^(0) - C xu (t)]. (43) 
7/ 

The two eigenvalues of the matrix A, determining the time-scales of the system 
1/A+ and 1/A_, read: 
1 



r ± ^/r 2 - Ar] . (44) 



As suggested by the interpretation given in Fig. El the parameters Tf = jf/k and r s 
should act as time-scales when they are separated enough. Indeed, an inspection of 
formula (1441) shows that the inverse of A + and A_ are proportional to 7/ and r s when 
they are well separated, i.e. r s ^> Tf or Tf ^> r s . However this is a limit case, and more 
general conditions can be considered. 

Our choices of parameters, always with Tf ^ T s , are resumed in Tabled! a case (a) 
where the time-scales are mixed, and three cases (b), (c) and (d) where scales are well 
separated. In particular, in cases (c) and (d), the position of the intermediate plateau is 
shifted at one of the extremes of the parametric plot, i.e. only one range of time-scales 
is visible. Of course we do not intend to exhaust all the possibilities of this rich model, 
but to offer a few examples which are interesting for the following question: what is the 
meaning of the usual "incomplete" parametric plot Xxx versus C xx , which neglects the 
contribution of Q xu 7 

The parametric plots, for the cases of Table (H are shown in Figure HI In Figure [51 
we present the corresponding contributions Q xx {t) and Q X u(t) as functions of time. We 
briefly discuss the four cases: 

(a) If the timescales are not separated, the general form of the parametric plot, see 
Fig. Hh, is a curve. In fact, as shown in Fig. [5h, the cross term Q xu {t) is relevant at 
all the time-scales. The slopes at the extremes of the parametric plot, which can 
be hard to measure in an experiment, are 1/Tf and Soo ^ 1/T S . Apart from that, 
the main information of the parametric plot is to point out the relevance of the 
coupling of x with the "hidden" variable u. 
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Figure 4. Parametric plots of integrated response Xxxit) versus self-correlation C xx (t) 
for the model in Eq. (|32[) with parameters given in Table [T] Lines with slopes equal to 
1/T S , l/Tf and Soo are also shown for reference. 



(b) In the "glassy" limit t s ^> Tf, with the constraint y = -f- ^ ~ 1/2, the well known 
broken line is found, see Fig. Hb, as discussed at the end of the previous section. 
Figure Eb shows that Q xu (t) is negligible during the first transient, up to the first 
plateau of x{^)i while it becomes relevant during the second rise of x{t) toward the 
final plateau. 

(c) If Tf ^> r s , the parametric plot, Fig. 0b, suggests an equilibrium-like behavior 
(similar to what one expects for Tf = T s ) with an effective temperature 1/soo 
which is different from both Tf and T s . Indeed, this case is quite interesting: the 
term Q xu is of the same order of Q xx during all relevant time-scales, but Q xu /Q xx 
appears to be almost constant. This leads to observe a KFR-like plot with a non- 
trivial slope. The close similarity between Q xx and Q xu is due to the high value of 
the coupling constant k' = -y s /r s . 

(d) In the last case, always with Tf 3> t s , the contribution of Q xu {t) is negligible at all 
relevant time-scales (see Fig. E]), giving place to a straigth parametric plot, shown 
in Fig. HI with slope 1/Tj. The low value of the coupling constant k' is in agreement 
with this observation. 

The lesson learnt from this brief study is that the shape of the parametric plot 
depends upon the timescales and the relative coupling k'/k. This is consistent with the 
fact that the correct formula for the response is always the GFR: ||^y = Q xx + Q xu - 

However, the definition of an effective temperature through the relation T e ff(t)^j^ = 
Q xx (t) in general (see case a), does not seem really useful. In particular limits, the 
behavior of the additional term Q xu is such that R oc Q xx in a range of time-scales, and 
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Figure 5. Integrated response Xxxif) as a function of time, for the model in Eq. (|32|) 
with parameters given in Table [T] The curves Q xx (t) and Q xu (t), representing the 
two contributions to the response, i.e. Xxx(t) = Qxx(t) + Qxu(t), are also shown. The 
violet curve with small circles represents the ratio Q X u(t)/Q xx (t)- 

therefore the measure of T e ff becomes meaningful. 
4. Granular gases 

In the previous Section we have shown in the analytically tractable case of linear 
Langevin equations how the presence of coupling among different degrees of freedom 
plays a role in the specific form of the GFR, which is in general different from the 
EFR and KFR. Such a feature is certainly not specific of the considered model: the 
non-equilibrium dynamics of a many particles system may present the same kind of 
nontrivial correlations among degrees of freedom, usually due to strong inhomogeneities 
generated by the lack of conservation laws valid at equilibrium. In these models the 
stationary distribution is not known and the use of the GFR for response analysis is 
rather subtle. 

In the following we shall analyze granular gases in the steady state, which offer an 
interesting benchmark of this ideas. 

4-1. The model 

Let us consider a <i-dimensional model for driven granular gases [27J [28j [29J E0] : N 
identical disks (in d — 2) or rods (in d — 1) of diameter 1 in a volume V = L x L (in 
d — 2) or total length L (in d — 1) with inelastic hard core interactions characterized 
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by an instantaneous velocity change 



v i = v ; y~ [( v i- v j 



■ n\n, 



(45) 



where i and j are the label of the colliding particles, v and v' are the velocity before 
and after the collision respectively, n is the unit vector joining the centers of particles 
and r e [0, 1] is the restitution coefficient which is equal to 1 in the elastic case. Each 
particle i is coupled to a "thermal bath" , such that its dynamics (between two successive 
collisions) obeys 



where r& and TJ, are parameters of the "bath" and rji(t) are independent normalized 
white noises. We restrict ourselves to the dilute or liquid-like regime, excluding more 
dense systems where the slowness of relaxation prevents clear measures and poses doubts 
about the stationarity of the regime and its ergodicity: in practice we consider packing 
fractions (fraction of occupied volume) if) = N/(4V) in the range 0.01 4- 0.5. Two 
important observables of the system are the mean free time between collisions, r c , and 
the so-called granular temperature T g = (|v| 2 )/<i. 

In this model is possible to recover two different regimes. When the thermostat is 
dominant, i.e. when a = Tc/t^ ^> 1, grains thermalize, on average, with the bath 
before experiencing a collision and the inelastic effects are negligible. This is an 
"equilibrium-like" regime, similar to the elastic case r = 1, where the granular gas 
is spatially homogeneous, the distribution of velocity is Maxwellian and T g = On 
the contrary, when r c < r&, non-equilibrium effects can emerge such as deviation from 
Maxwell-Boltzmann statistics, spatial inhomogeneities and T g < T b [271 [28j [29j [30] . This 
"granular regime", easily reached when packing fraction is increased or inelasticity is 
reduced, is characterized by correlations among different particles. This peculiarity is 
the key ingredient for a correct response analysis of these systems, as we shall see in the 
following. 

4-2. Failure of the EFR for strong dissipation 

An analysis of the FDR for the previous model has been performed in [311 [32] (in d = 2) 
and [33J (in d = 1), and discussed also in [9]. Similar results are also obtained for other 
models, such as the inelastic Maxwell model on a d = 2 lattice driven by a Gaussian 
thermostat and mean field granular gases [3U [35]. The protocol used in numerical 
experiments cited above is the following: 

(i) the gas is prepared in a "thermal" state, with random velocity components extracted 
from a Gaussian with zero average and given variance, and positions of the particles 
chosen uniformly random in the box, avoiding overlapping configurations. 

(ii) The system is let evolve until a statistically stationary state is reached, which is 
set as time 0: we verify that the total kinetic energy fluctuates around an average 
value which does not depend on initial conditions. 




(46) 
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(iii) A copy of the system is obtained, identical to the original but for one particle, 
whose x (for instance) velocity component is incremented of a fixed amount 5v (0). 

(iv) Both systems are let evolve with the unperturbed dynamics. For the random 
thermostats, the same noise realization is used. The perturbed tracer has velocity 
v'(t), while the unperturbed one has velocity v(t), so that 5v(t) = v'{t) — v(t). 

(v) After a time t max large enough to have lost memory of the configuration at time 0, 
a new copy is done with perturbing a new random particle and the new response is 
measured. This procedure is repeated until a sufficient collection of data is obtained 

(vi) Finally the autocorrelation function C vv (t) = (v(t)v(0)) in the original system and 
the response R vv (t) = are measured. 

The main result of those studies can be so summarized: in the dilute limit (where 
the packing fraction ip — > 0) or in the elastic limit (r — > 1), or in the limit of efficient 
thermostat a — >• oo (which is usually implied by the dilute limit), the Einstein relation, 
Eq. ([9]), is recovered for the velocity of a tracer particle, i.e. R vv (t) = C vv (t)/T g . On 
the contrary, when these conditions fails one can observe strong deviations from EFR. 
In addiction, there are same remarkable points: 

• Non-Gaussian velocity distributions can appear also in the dilute regime, but they 
seems to have a minor role in the violations of EFR [32J . 

• The same scenario holds in dimension d — 1, where the tracer is sub-diffusive, i.e. 
where (\x(t) — x(Q)\ 2 ) ~ t 1 ! 2 which also implies a non-monotonic C vv {t) with a 
power-law tail C vv (t) ~ — 1~ 3 / 2 for large t [33J. 

• When a mixture of two different kinds of grains (e.g. with different masses 
or different restitution coefficients) is considered [36J, the two components bear 
different granular temperatures and this lead to separated FDR in the dilute limit, 
i.e. a tracer satisfies the EFR with its own temperature, making very difficult to 
obtain a neutral thermometer based on FDR. 

The violation of the EFR is more and more pronounced as the inelasticity increases 
(lower values of r), the importance of the bath is reduced (lower values of a) or the 
packing fraction ip is increased, as shown in Figure [6] In correspondence of such 
variations of parameters, the correlation between velocities of adjacent particles is also 
enhanced, a phenomena which is ruled out in equilibrium fluids. This effect can be 
directly measured in many ways, a possibility is shown in Figure [7J In conclusion, the 
general lesson is that there is a quite clear correspondence between violations of the 
EFR and the appearance of correlations among different degrees of freedom. 

4-3. A correct prediction of the response 

Let us begin with an example explaining how a blind comparison between 
autocorrelation and response can be misleading. Clearly the EFR is satisfied in the 
elastic ("equilibrium") case, i.e. R vv (t) oc C vv (t), thanks to the fact that, in the invariant 
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Figure 6. Parametric plots to check the EFR, for d = 1 and d = 2 models of inelastic 
hard-core gases with thermal bath. Different choices of parameters r (restitution 
coefficient), a = T c /rb and ip (packing fraction) are shown: note that one can change 
a at ip or r fixed (changing Tb), but - in general - changes in tp or r determine also 
changes in a (because of changes in r c ). In all plots, the dashed line marks the Einstein 
relation R vv — C vv (t) / C vv (0) . 
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Figure 7. Static (same time) particle-particle correlations for energy fluctuations, 
for the d = 1 inelastic hard rods gas. In the inelastic case (squares) the coefficient 
is higher when the dissipation gets stronger, i.e. for small a (left) or high packing 
fraction (right). The coefficient in the elastic case is also shown, which is negligible for 
all the values of the parameters. 
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Figure 8. Comparison of an clastic (left) and inelastic (right) gas of hard disks 
coupled to the thermal bath in d = 2: the response R vv (denoted as R for simplicity), 
the susceptibility Xwit) = J R V v(t')dt' and the normalized velocity auto-correlation 
C(t) — C vv (t)/C vv (0) are plotted in several different ways. The top graphs show the 
validity or breakdown of the Einstein relation R — C for the elastic or inelastic case 
respectively; the bottom graphs display the parametric plot Xw vs. C vv which would 
follow the form Xw = 1 — C if the second form of FDR, Eq. (fTTj) , i.e. R vv = —C vv , 
held. 



measure, there is a very weak coupling between the tracer velocity and other degrees 
of freedom. However, it is also interesting to note that the shape of C vv (t) is far from 
being an exponential, also in the elastic ( "equilibrium" ) case, because of the presence 
of two characteristic times rj and r c . This non-exponential shape of C vv (t) leads to a 
failure of the KFR, Eq. ([12]) . as put in evidence in the lower left plot of Fig. [HI It is also 
instructive to plot the same quantities for a strongly inelastic case, see right plots of 
Fig. [HI where an almost linear relation seems at work. A first rough explanation is the 
lower granular temperature which is responsible for a higher mean free time r c , implying 
that t c and tj, are slightly closer to each other with respect to the elastic case, making 
the C vv (t) similar to a single exponential. 

This example show how the only way to have a correct prediction of the response 
can reside in the use of the GFR, which is always valid. For a quantitative comparison 
between correlation functions and response, one needs some hypothesis on the stationary 
probability distribution, in particular about the kind of coupling between different phase- 
space variables. We report the result of a simple assumption, partly inspired to an idea 
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of Speck and Seifert [37], where correlation among variables is mediated by a fluctuating 
"hydrodynamic" velocity field u(x), in such a way that the relevant part of the stationary 
probability distribution for the tracer reads, approximately: 

P m (v,x, f )~exp{- ' V - 2 "^ f >l 2 }, (47) 

with u(x, t) a local velocity average, defined on a small cell of diameter L^ ox centered in 
the particle. This is motivated by the observation that, at high density or inelasticities, 
spatially structured velocity fluctuations appear in the system for some time, even in the 
presence of external noise [201 EE]- The generalized FDR following from assumption (T4TI) 
reads 

R m (t) = C.(t) = ±-{v{t){v(0) - u[x(0)}}). (48) 

and is nicely verified in Fig. [9j In other words, one has a correction to the "naive" 
expectation R vv (t) = (v(t)v (0))/T g , i.e. the extra term —(v(t)u[x(0)]) originated by 
the presence of a "hydrodynamic" velocity field. A recent experiment [39] shows the 
presence, in a rather clear way, of a similar extra term (with respect to the KFR) for a 
colloidal particle in a toroidal optical trap, in a non-equilibrium steady state. 
We conclude this section underlining the connections between the results obtained for 
granular gases and Langevin equations. In general, the behaviour of the Langevin model 
and of the granular model show some differences, for example the case expressed in 
Figure 0] (case b) has not a counterpart in these models and the "effective temperature" 
approach is meaningless, even when times scales are well separated. However the use 
of GFR shows how, in both examples, the response is given by a sum of different 
contributes, and in some special limits, the cross correlation term can be neglected 
and a comparison between the response and the autocorrelation does make sense. This 
happens in the "equilibrium-like" cases of the Langevin model (cfr. Figure HI cases c,d) 
and in the dilute regime of the granular gas, in which there is no coupling between the 
velocity of the tracer and the "hidden variable" embodied by the local velocity average. 
On the contrary, when this approximation is not correct, a response-autocorrelation plot 
shows strong deviations from linearity, but can be predicted by taking into account all 
the contributes in the computation of the response. 

5. Conclusions 

In this paper we have reviewed different forms of the fluctuation-dissipation relation for 
steady states. The most general one is the GFR, Eq. d^J), which requires the knowledge 
of the relevant degrees of freedom and their reciprocal couplings in the system. When 
this knowledge is not accessible, the study of the response to a perturbation of a certain 
variable, compared to the correlation of that variable in the unperturbed state, has not 
a simple meaning, in general. 

As an example, we study in detail two limits of a generalized Langevin equation 
with memory, for a particle which moves in a harmonic potential and is in contact 
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Figure 9. Response and correlations of the velocity tracer in the d = 2 gas of 
inelastic hard disks for a choice of parameters such that FDR type (i) (Einstein 
relation, Eq. ((9])) is strongly violated. The dashed green curves show the conjecture 
R vv = C s {t) = Jh(v(i){v(0) — u[x(0)]}} is probed, where u(x) is the local velocity field 
measured by coarse-graining in boxes (centered with the tracer) of size Lbox- 

with two thermostats at different temperatures. In the overdamped case, the response- 
correlation parametric plot may reveal a broken line shape, where the two slopes are 
given by the inverse temperatures of the baths: a necessary, but not sufficient, condition 
for this to happen is that time-scales are well separated. In the general case, however, 
the plot can be more difficult to read, showing intermediate "effective temperatures", 
as well as a more general non-linear shape. 

Nevertheless, the problem can be recast as a Markovian dynamics, by means of 
the introduction of additional degrees of freedoms. The Markovian dynamics has a 
stationary distribution and satisfies a general relation (GFR) between the response and 
a specific correlation function, which is related to the stationary measure and takes 
into account the coupling between degrees of freedom. We show how this coupling is 
responsible for the "violation" of the usual FDR. 

An interesting case where correlations play a role in the "violation" of FDR, 
is the dynamics of a tracer particle in a driven granular gas. Here the response 
is not proportional, in general, to the unperturbed autocorrelation and an effective 
temperature does not seem to be informative, even when the time-scales are well 
separated. Nevertheless a GFR should always be valid, provided an appropriate 
description of the dynamics is given. Indeed, as already stressed by Onsager and 
Machlup in their seminal work on fluctuations and irreversible processes [40] . a basic 
ingredient for a "good statistical" description, is Markovianity. We remind their 
important caveat: how do you know you have taken enough variables, for it to be 
Markovian? We have seen that, following this suggestion, it is possible to have deeper 
insight on the so-called "violations of FDR" and in some cases (e.g. in Fig. [9]), one can 
try to guess the correct correlation involved in the dynamics. 
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